Revisit the ACDM Universe in f(R) gravity 
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bJQi Abstract 

We perform a number of explicit reconstructions for the f(R) gravity model to give rise to the 
particular ACDM evolution of the universe. We find well-defined real valued analytical forms for 
the f(R) model to describe the universe both in the early epoch from radiation to matter dominated 
era and the late time acceleration period. We further examine the viability of the derived f(R) 
Q-i model and find that it is viable to describe the universe evolution in the past and there does not 

exist the future singularity in the Lagrange. 
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There has been accumulated conclusive evidence from supernovae[l] and other 
observations [2, 3] in the last decade, indicating that our Universe is undergoing a phase 
of accelerated expansion. Understanding the nature of the cosmic acceleration is one of the 
biggest questions in modern physics. The leading explanation of the accelerated expansion 
is the cosmological constant in the context of General Relativity. However, the measured 
value of cosmological constant is far below the prediction of any sensible quantum field theo- 
ries and furthermore the cosmological constant leads inevitably to the coincidence problem, 
namely why the energy densities of matter and the vacuum are in the same order today(see 
[ I] for review). 

Alternatively, another possible explanation of the acceleration is the modified gravity. 
The General relativity might not be ultimately correct in the cosmological scale. The Uni- 
verse might be described by some kind of modified gravity. One simplest attempt is called 
f{R) gravity, in which the scalar curvature in the Lagrange density of Einstein's gravity is 
replaced by an arbitrary function of R. However the complexity of the field equations makes 
it difficult to obtain a viable f(R) model to satisfy both cosmological and local gravity 
constraints [5]. Recently, there appeared a useful approach to reconstruct the f(R) model 
from the observed expansion history of the universe and try to invert the field equations 
to deduce what class of f(R) theories give rise to the particular cosmological evolution 
[6] [7] [8] [9] [10] [11]. Some analytical forms for f(R) gravity to admit the ACDM expansion 
history in the background spacetime were constructed [9] [10]. However, it was argued that 
only simple real expression of f(R) model in the Lagrange could admit an exact ACDM 
expansion history [10]. 

In this paper, we will further study this problem. We will perform a number of explicit 
reconstructions which lead to a number of interesting results. We will show that we can 
derive a well-defined real analytical f(R) form in terms of the hypergeometric functions to 
admit an exact ACDM expansion history. We will disclose that the f(R) gravity not only 
can admit an exact ACDM expansion in the recent epoch of the Universe but also can admit 
an exact ACDM expansion in the early time of the Universe. 

We consider the background dynamics of the Universe in f(R) gravity described by[12] 

H* = ^—±-H^ + ^p . (1) 
6F F 3F H v ; 

where p = p m + p r , F = gjj/ and k 2 = 8ttG. If we convert the derivative from the cosmic 



time t to x = In a and further take the derivative of the above equation, we obtain 

+ f l dlnE + ^ f = ( 2 ) 

dx 2 2 dx dx dx 3E dx ' 



where 



* S 3(| + 4E) , (3) 

For convenience, we take the energy density p and the scalar curvature R in Eq.(2) in the 
unit of Hq and we also set k 2 = 1 in our analysis. In order to get a viable f(R) model with a 
reasonable expansion history of the Universe allowed by observations, we can parameterize 
E(x) in Eq.(2) as the standard model in Einstein's gravity with an effective dark energy 
equation of state(EoS) w[6][7] 

E(x) = Q°e- 4x + Q° m e- 3x + Q^s f^i+u,)dx ^ ^ 
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After specifying the expansion history of the Universe E(x), Eq.(2) becomes a second 
order differential equation of E(x). If we can find the solution of Eq.(2), we then obtain 
the explicit form of f(R) correspondingly We find that it is more convenient to use the 
quantity G = F — 1 instead of F, so that Eq.(2) can be changed into 

d*G | ldlnJg dG | d\nE c _ 3(l + w)n° d c _ 3f * a+w)dx _ (g) 
dx 2 2 dx dx dx E 

To mimic the exact ACDM expansion history of the universe with w = —1, Eq.( 6) is a 
homogenous equation. We will focus on this case hereafter and show that the solution of 
the above differential equation will directly lead to the real valued f(R) form which gives 
rise to the cosmological evolution as that of the ACDM model. 



Eq.(6) does not have the overall analytical solution with the full expression of Eq.(4) 
to describe the universe from the radiation dominated epoch to the late time acceleration. 
However, Eq.(6) does have analytical solutions in different epochs in the evolution of the 
Universe. Let's first concentrate on the early evolution of the Universe from the radiation 
dominated epoch to the matter dominated epoch. In this case, E can be taken as 

E x ~ n° m e- 3x + n° r e~ 4x , 

and the general solution of Eq.(6) has the form 

(U- n± 0° \ 0° 
l"7^ eX ' )+D 1 e Ax 2 F 1 [m.,m + ;6;~-^e x } (7) 
— 1 4 r / r 

where G22 is the Meijer G function, 2 -Pi is the Gaussian hypergeometric function and C±, 
D\ are arbitrary constants which can be determined by boundary conditions. 
The indexes in the solutions are 
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The viable f(R) models should be in a "chameleon" type [13] [14] which provides a mechanism 
to pass the local test. However, the first term of Eq.(7) is divergent when x goes to the 
minus infinity. Thus the requirement lmx E _ i ._ 00 G(x) = puts the condition C\ = 0, so that 
G\{x) turns out to be 

Gt(x) = D 1 e Ax 2 F 1 [m.,m + ;Q;-^e x }. (8) 

We can obtain the explicit form for the f(R) model to describe the early universe by doing 
the integration 

G(x)-^dx , (9) 
where the scalar curvature R can be written as 

i?^3^e- 3 * , (10) 
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which only contains the component of matter because the radiation does not have any 
contribution to the scalar curvature R since its energy momentum tensor is traceless. Eq.(lO) 
is valid in the expansion history of the Universe from the radiation dominated epoch to the 
deep matter dominated epoch. The term 4^ in Eq.(9) thus can be expressed as 



dR 

dx 



-^° m e- 3x 

where x, in turn, can be presented in terms of R 

x(R) = 



When x goes to the minus infinity x — > — oo, R goes to infinity R 
above equations, we obtain the explicit expression for f(R) as 



(11) 
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(12) 

where fi(R) and R are in the unite of Hq. Eq.(12) is mathematically well-defined for 
all the positive value of the scalar curvature R > because the hypergeometric function 
2-F\ [a, b; c; z] has the integral representation on the real axis when b > and c > 



jFi[a,6; c; z] 



V{b)Y{c -b) X J 1 



l-t) c -"- l (l - zt)~ a dt 



(13) 



where T is the Euler Gamma function. The above expression is well-defined in the range 
— oo < z < 1 and the resulting value of 2F\[a, b; c; z\ is also a real value. In order to present 
the f(R) in SI units, we can insert Eq.(5) and then we obtain 



fi(R) = R 
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(14) 



where p° and are the energy density of radiation and matter at today respectively. Rq is 
the scalar curvature Rq = K 2 Pm + K Vr + 4A. Eq.(14) shows that we have the well-defined 
real analytical function for f(R) gravity that can exactly reproduce the same background 
expansion as that of the ACDM model from the radiation dominated epoch to the matter 
dominated epoch. At very early time when the curvature is very high, R 3> 4A, we have 
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Thus, Eq.(14) reduces to 



(K 2 fl 0\ V3 



R 

-2^0 



(16) 



which shows that in the high curvature region R » u pj,, fi(R) goes back to the standard 
Einstein's model fi(R) ~ R. 

Next we turn to investigate the most interesting case and try to answer whether the f(R) 
model can mimic the evolution of the Universe from the matter dominated epoch to the late 
time acceleration. In this case, E can be taken as 

e 2 = n° m e- 3 * + n° d , 

where is a constant. The scalar curvature can be presented as 

R = 3fi° n e" 3x + 1202 • 
The general solution of Eq.(6) gives 
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where the indexes are 
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After doing integration, we can get the explicit expression for f(R) from Eq.(9) 

f 2 (R) = R-6Q° d + e.+e + , 

where 
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Inserting Eq.(5), we obtain 
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where w x = D 2 (R - AA)p+/( P+ - 1)/A P+_1 and ro 2 = C 2 (R - AA) p -/(p_ - l)/A p -~\ The 
constant parameter A is defined as A = n 2 pd and pa is the effective energy density of dark 
energy. When w\ = w 2 = 0, A is just the cosmological constant. Noting the fact that 
p_ — 1 < 0, when R — > +00, the second term of Eq.(20) becomes divergent. This is not 
allowed for the "chameleon" type solution, so that we need to set w 2 = 0. Therefore the 
solution has the form 



f 2 (R) = R-2A-w l 
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From the integral representation of the hypergeometric function Eq.(13), it is clear that 
Eq.(21) is mathematically well-defined in the range R > 4A. Eq.(21) is a real function in 
the physical range from the matter dominated epoch to the future expansion of the Universe. 

Eq.(6) is our starting point to find the analytic expression for the f(R) model to mimic 
the ACDM cosmology. From Eq.(17), we can see clearly that the scalar curvature R obtained 
from Eq.(6) can be constrained automatically in the physical range R > 4A. Our starting 
point is different from that in [9] [10], where they got their solution by solving the differential 
equation [9] [10] [11] 



- 3(J? - 3A)(R - 4A)^/(fl) + (| - 3A)^/(i2) + \f{R) +4A-R 







(22) 



This equation was derived from Eq.(l). In Eq.(22), R can be chosen as any value on the 
real axes, so that not all solutions of Eq. (22) are physical. We need carefully to analyze the 
solutions of Eq.(22) . 

A particular solution for Eq.(22) is 



f p (R) = R-2A 



(23) 



Therefore, we only focus on the homogenous solutions of Eq.(22) hereafter since f(R) = 
f p (R) — fh(R)- The homogenous part of Eq.(22) is a standard hypergeometric equation 



which may have at most 24 solutions in the complex plane around three different singular 
points(i? = oo,3A,4A) [15]. However, if we focus on real solutions, Eq.(22) may have at 
most 32 solutions on the real axes around four different singular points R = — oo, 3A, 4A, +oo 
. We will extensively discuss all of these solutions in the following. 
The solution around +oo reads, 
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(24) 



The above expression is just Eq.(20) which is the physical solution of Eq.(22). We can see 
clearly that when R — > +oo, fh{R) is well-defined on the real axes and actually fh{R) is a 
real function for the whole range of R > AA as discussed previously. R = AA is a finite point 
because lim R ^ 4 A fh(R) is a finite value. However, the derivatives jj}fh(R) of all the terms 
in the above expression are divergent at R — A A. 
The solution around — oo reads 
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(25) 



This expression was obtained in [9]. This solution is well-defined when R — > —oo. It is a 
real function when R < 3A. However, it becomes complex when R > 3A. R = 3A is a finite 
point. fh(R) and its derivative jj^fh(R) are well defined at R = 3A. Clearly Eq.(25) is not 
a physical solution, it can not satisfy Eq.(6). 
The solution around 3A reads, 



f h (R) = w 12 F, 



1 R 
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where a± = (—7 ± v73)/12 and /3± = (11 ± \/73)/12. The above solution was derived in 
[10]. fh(R) and its derivative jjifh(R) are well-defined on the finite singular point R = 3A. 
When R > AA or R < 3A, the second term of Eq.(26) becomes complex. However, in 
contrast to what was claimed in [10], when R > AA the homogenous part of Eq.(22) do have 
the real analytical solution, namely Eq.(24). 



8 



The solution around 4A reads, 
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Clearly the above solution is well-defined on the finite singular point R = 4A and the solution 
is valid for R > 4A. It is a solution in the physical range and satisfies Eq.(6). Eq.(27) is 
equivalent to Eq.(24) since they are defined in the same range. Eq.(27) is simply a new linear 
combination of the solutions in Eq.(24). However, Eq.(27) has different behaviors at singular 
points if compared with Eq.(24) owing to the different linear combination of hypergeometric 
functions. For instance, the derivatives j^fh(R) of all the terms in Eq.(24) are divergent at 
R = 4A while in Eq.(27), the derivative of the first term 2F1 \p + — 1, j>_ — 1; |; — — 4)] 
is well-defined at R = 4A. Furthermore, Eq.(27) apparently does not have the "chameleon" 
property because both terms in Eq.(27) are divergent when R goes to infinity. However, their 
linear combination, namely, the term in Eq.(21) is convergent when R — > +00. Therefore, 
even for equivalent results, we should carefully choose the proper expressions according to 
boundary conditions at singular points. 

Using the Euler transformation and Pfaff transformation, from Eqs. (24,25,26,27), we 
can find out all the 8 x 4 = 32 real solutions for Eq.(22). Eqs. (24,25,26,27) complete 
the different behaviors at different singular points for the solutions of Eq.(22). Although 
Eqs. (24,25,26,27) are substantially different on the real axes, when extended to the com- 
plex plane, Eqs. (24,25,26,27) are equivalent to each other because they can be related by 
connection formulas. However, for different expressions, they have different behaviors at 
the singular points. It should be very careful to choose the proper expressions accord- 
ing to different boundary conditions. The physical solutions for f(R) models to describe 
the Universe should be well-defined in the range R > 4A and possess the "chameleon" 
property \im R _> +OQ f h (R) = 0, \im R ^ +OQ j^f h (R) = 0. The only physical solution is the 
"chameleon" part of Eq.(24)(o7 2 = 0), namely Eq.(21). The result obtained from Eq.(22) 
has bigger range in the solutions than that obtained from Eq.(6). We need very carefully 
to pick up the physical solution. It is more convenient to start from Eq.(6) to find physical 
solutions. 

Having the well-defined analytical expression for f(R) model, we will further discuss its 
viability. In the very early Universe, fi(R) ~ .M-R) ~ R- In order to evade the instabilities 
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Figure 1: From top D 2 = 0.05, 0.03, 0.01, 0, -0.01, -0.03, -0.05 
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Figure 2: From top D 2 = -0.05, -0.03, -0.01, 0, 0.01, 0.03, 0.05 

of f(R) model, it requires that F > and f RR > [12] [16]. From Fig.(l) and Fig.(2), 
we can see that when D 2 < 0, in the past expansion of the Universe x < 0, we have 
F 2 (R) > 0,f 2 Rn{R) > so that the model f 2 {R) is viable in the past. However, for 
the future expansion of the Universe, the condition D 2 < can not guarantee F 2 (x) to 
be always positive. The zero-crossing behavior will lead to singularities in the conformal 
transformation [ J 7] and the negative F will lead to the imaginary mass of particles m = 
m/y/F in the Einstein frame [12] [18]. Furthermore, we can see from Fig.(l) and Fig.(2) 
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that, when x — > +00 or R — > 4A, the derivative of f(R), namely, F(R) = j^f(R) and 



fRR,(R) — dR~F(R) are divergent. Although the expression of f2(R) has such weakness 
in the future evolution of the universe, /^(-R) does not have the future singularity in the 
Lagrange because /^(-R) is finite at R = 4A. The future point happens at 



x— >+oo 



lim i? = 4A , (28) 

X ' 

and f2(R) is finite at R = 4A 



lim M R) = 2A - w 4^511 + 79V73)r(2/3)r(-r„) « 2A - 1.256* , 

^ (-5 + v / 73)(-l + v / 73)(7+v / 73)r(-j9_)r(g + ) 

(29) 

when w < 0, we can find that ^(4A) > 2A. 

In summary, in this work we have constructed the f(R) gravity model to mimic the 
ACDM universe expansion in the early and late epochs. We found that there exists a real 
valued function of Ricci scalar presented in terms of hypergeometric function, which can give 
rise to the particular cosmological evolution as the ACDM model. Although the constructed 
f{R) model has weakness in describing the future expansion of the Universe in the Einstein 
frame, in the Jordan frame it is viable to describe the past evolution of the Universe and it 
does not have the future singularity in the Lagrange. 



Acknowledgments 

J.H.He acknowledges the Financial support of MIUR through PRIN 2008 and ASI through 
contract Euclid-NIS 1/039/10/0. The work of B.Wang was partially supported by NNSF 
of China under grant 10878001 and the National Basic Research Program of China under 
grant 2010CB833000. 



[1] S. J. Perlmutter et al., Nature 391 (1998) 51; A. G. Riess et al., Astron. J. 116 (1998) 1009 ; 
S. J. Perlmutter et al., Astroph. J. 517 (1999) 565 ; J. L. Tonry et al., Astroph. J. 594 (2003) 
1; A. G. Riess et al., Astroph. J. 607 (2005) 665 ; P. Astier et al., Astron. Astroph. 447 (2006) 
31 ; A G. Riess et al, Astroph. J. 659 (2007) 98. 

[2] E. Komatsu, et.al., ApJS 192 (2011) 18, arXiv:1001.4538. 



11 



[3] A. G. Sanchez, et.al. arXiv:1203.6616. 

[4] S. M. Carroll, LivingRev. Rel. 4 (2001) 1, astro-ph/0004075. 

[5] P.G. Bergmann, Int. J. Theor. Phys. 1 (1968) 25;A. A. Starobinsky, Phys. Lett. B 91 (1980) 
99; A.L. Erickcek, T.L. Smith, T.L., and M. Kamionkowski, Phys. Rev. D 74 (2006) 121501; 
V. Faraoni, Phys. Rev. D 74 (2006) 023529; S. Capozziello, and S. Tsujikawa, Phys. Rev. D 
77 (2008) 107501 ;T. Chiba, T.L. Smith and A. L. Erickcek, Phys. Rev. D 75 (2007) 124014; I. 
Navarro, and K. Van Acoleyen, J. Cosmol. Astropart. Phys. 02 (2007) 022; G. J. Olmo, Phys. 
Rev. Lett. 95 (2005) 261102;G. J. Olmo, Phys. Rev. D 72 (2005) 083505; Amendola, L., Po- 
larski, D., and Tsujikawa, S., Phys. Rev. Lett. 98 (2007) 131302, astro-ph/0603703; L. Amen- 
dola, R. Gannouji, D. Polarski, S. Tsujikawa, Phys.Rev.D 75 (2007) 083504, gr-qc/0612180; 
L. Amendola, Phys.Rev.D 60 (1999) 043501, astro-ph/9904120; S. Nojiri, Sergei D. Odintsov, 
Phys.Rev.D68 (2003) 123512; S. Nojiri, S.D. Odintsov, Int.J.Geom.Meth.Mod.Phys. 4 (2007) 
146; 

[6] Y-S Song, W. Hu, and I. Sawicki, Phys.Rev. D 75 (2007) 044004 ,astro-ph/0610532; Y-S Song, 

H. Peiris, W. Hu, Phys.Rev.D76 (2007) 063517,arXiv:0706.2399. 
[7] L. Pogosian, A. Silvestri, Phys.Rev.D77 (2008) 023503, arXiv:0709.0296. 
[8] L. Lombriser, A. Slosar, U. Seljak, W. Hu, arXiv: 1003.3009. 

[9] A. de la Cruz-Dombriz, A. Dobado. Phys.Rev. D74 (2006) 087501, arXiv:gr-qc/0607118. 
[10] P. K. S. Dunsby et al., Phys.Rev. D82 (2010) 023519, arXiv: 1005.2205. 

[11] T.Multamaki and I. Vilja, Phys. Rev. D73 (2006) 024018,astro-ph/0506692; S. Nojiri and 

S. D. Odintsov, Phys. Rev. D74 (2006) 086005, hep-th/0608008; Shin'ichi Nojiri, Sergei D. 

Odintsov J.Phys.A40 (2007) 6725, hep-th/0610164; S. Capozziello, S. Nojiri, S. D. Odintsov, 

A. Troisi, Phys. Lett. B639 (2006) 135, astro-ph/0604431. 
[12] A. Silvestri, M. Trodden, Rept.Prog.Phys.72 (2009) 096901, arXiv:0904.0024; A. Felice, S. 

Tsujikawa, LivingRev. Rel. 13 (2010) 3 , arXiv: 1002.4928; T. Clifton, P. G. Ferreira, A. 

Padilla, C. Skordis, Physics Reports 513 (2012) 1 ; T. P. Sotiriou, V. Faraoni, Rev. Mod. 

Phys. 82 (2010) 451, arXiv:0805.1726; S. Nojiri, Sergei D. Odintsov, Phys.Rept.505 (2011) 59, 

arXiv: 101 1.0544. 

[13] D. F. Mota, J. D. Barrow, Phys.Lett.B581 (2004) 141, astro-ph/0306047. 
[14] J. Khoury, and A. Weltman, Phys. Rev. D 69 (2004) 044026 ; J. Khoury, and A. Weltman, 
Phys. Rev. Lett. 93 (2004) 171104. 



12 



[15] M.Abramowitz, I.Stegun, Handbook of Mathematical Functions, Dover (1970). 

[16] I. Sawicki, W. Hu, Phys. Rev. D 75 (2007) 127502, arXiv:astro-ph/0702278. 

[17] K. i. Maeda, Phys.Rev.D 39 (1989) 3159. 

[18] Jian-Hua He, Bin Wang, E. Abdalla, Phys. Rev. D 84 (2011) 123526, arXiv:1109.1730. 



13 



